Design of Low Band Gap Double Perovskites from First Principles 
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Using density functional theory (DFT)-based calculations, we propose a family of metastable, as- 
yet unmade V 5+ and Cr 6+ double perovskite compounds with low band gaps spanning much of the 
visible region of the solar spectrum. Through analysis of a related set of measured optical gaps of 
d° ABO3 perovskites and A2B'B06 double perovskites, an ad hoc procedure is developed to correct 
DFT and many-body perturbation theory gaps, bringing them into quantitative agreement with 
experiment for measured compounds, and predicting that V 5+ and Cr 6+ double perovskites would 
have gaps ranging from approximately 1.1-2.4 eV, significantly lower than previous materials studied 
in this class. DFT calculations also establish that these V 5+ and Cr 6+ compounds are likely able 
to be synthesized, either in bulk form or as epitaxial thin films. These compounds would comprise 
a new class of semiconducting double perovskites for potential use in solar energy conversion and 
other optoelectronic applications. 
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INTRODUCTION 

The ability to tunc the electronic structure of per- 
ovskite oxides is of significant interest in areas ranging 
from superconductivity to thermoelectricity to solar en- 
ergy conversion. For example, perovskite-derived transi- 
tion metal oxide compounds have been shown to photo- 
catalyze the splitting of water into hydrogen and oxygen 
gas. [ll, [H However, their ability to harvest sunlight in 
the visible spectrum is primarily limited by their large 
optical band gaps. Experimental efforts to reduce per- 
ovskite optical gap s have included chemical doping and 
substitution. [3l-llC| Especially relevant to this work, Eng 
et al. [lH synthesized a number of d° perovskites (ABO3) 
and double perovskites (A 2 B'B06) with rock-salt B-site 
ordering. Still, none of these compounds possess gaps 
below 2.7 eV, leaving open the question of whether there 
might exist materials in this class with more significant 
overlap with the solar spectrum. 

First-principles calculations based on density func- 
tional theory (DFT) have proven useful for exploring the 
large structural parameter space of chemically doped and 
substituted perovskites to search for new materials with 
reduced band gaps. 0, 0, Il2l - fl5| However, since DFT is 
well known to underestimate band gaps, the quantitative 
nature of these predictions is limited. Many-body per- 
turbation theory within the GW approximation can re- 
sult in significantly more reliable gaps for many inorganic 
semiconductors, [ly, [l7| despite its neglect of excitonic 
effects and electron-phonon interactions. Yet standard 
implementations of the GW approximation — in particu- 
lar one-shot GqWq corrections — have been less successful 
for transition metal oxides where d states feature promi- 
nently. (l7l - [2"5| Thus, quantitative first- principles predic- 
tion of band gaps for transition metal oxides remains a 
principal challenge. 

In this work, we propose a class of new double per- 



ovskite transition metal oxides with band gaps well 
within the visible region of the solar spectrum, mitigat- 
ing difficulties in using conventional DFT-based methods 
to quantitatively assess transition metal oxide band gaps 
with an empirical approach. We first identify a group 
of as-yct unmade d° double perovskites containing V 5+ 
or Cr 6+ at the B sites that are metastable with respect 
to structural distortions and changes in B-sitc ordering. 
Building on existing optical data on gaps of 30 differ- 
ent d° perovskites and double perovskites, [ll[ we then 
develop and use an ad hoc procedure, based on a fit to 
the difference between DFT and experimental gaps that 
corrects for DFT errors in electron correlation and ex- 
change, to predict optical gaps to within about ±0.2 eV. 
These ad hoc corrections are used to predict optical gaps 
of metastable V 5+ and Cr 6+ compounds in the middle 
or toward the lower end of the visible region of the so- 
lar spectrum, establishing these new double perovskites' 
potential for solar energy conversion and other optoelec- 
tronic applications. 



COMPUTATIONAL DETAILS 

All crystal structures and total energies are com- 
puted using DFT within the local density approxima- 
tion (LP A ), t he VASP package, [2^ and PAW poten- 



tials [27J. [55j A plane- wave cutoff of 500 eV is used 



throughout. Calculations on 5-atom perovskite unit cells 
use a T-centered 10 x 10 x 10 fc-point mesh, while larger 
unit cells use correspondingly fewer fc-points. Zone- 
center phonon frequencies are computed using density 
functional perturbation theory. Elastic constants and 
phonon dispersion relations are calculated using a finite 
differences approach and the Phonopy code. [28[ GqWq 
corrections to DFT-LDA eigenvalues (referred to here- 
after as GqWo@LDA) are evaluated in the LDA opti- 
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mizcd crystal structure, using a frequency-dependent di- 
electric function calculated within the random phase ap- 
proximation, sampling 100 frequency points; 144 total 
bands per 5-atom perovskite unit cell (approximately 20 
of which are occupied, depending on chemical composi- 
tion); and a plane- wave cutoff for the response function 
of 200 eV. As an additional point of comparison, hybrid 
HSE06 band gaps are calculated using HSE06 opti- 
mized crystal structures. Both G W @LDA and HSE06 
calculations on 5-atom perovskite unit cells use a T- 
centered 6x6x6 fc-point mesh, while larger unit cells 
use correspondingly fewer fc-points. 



RESULTS AND DISCUSSION 



leaving open the question of how much these chemical 
substitutions might reduce the optical gap relative to 
SrTiOa. A first estimate of their optical gaps (Fig. [Th) is 
quite promising: assuming that their LDA, GoWo@LDA, 
and HSE06 gaps have the same errors as for SrTiOa rela- 
tive to room-temperature absorption experiments (i.e., 
-1.30, 0.76, and 0.21 eV for LDA, G W @LDA, and 
HSE06, respectively), all twelve compounds are expected 
to have significantly smaller gaps than SrTiOa that span 
considerable regions of the solar spectrum, ranging from 
2.8 eV down to 0.8 eV in these initial estimates. How- 
ever, these estimates are severely limited, particularly 
given that LDA, HSE06, and GW are not rigorous theo- 
ries of room-temperature optical gaps; additionally, they 
carry unknown uncertainties. 



Selection of novel perovskite-derived compounds 

Based on measurements of 31 d° perovskite-derived 
compounds, Eng et al. ll| rationalized that their op- 



tical band gaps depend primarily on the "effective elec- 
tronegativity" of the B-site transition metal cation. That 
is, smaller, more highly charged metal cations result 
in lower-energy conduction bands and correspondingly 
smaller gaps. We use this principle to infer that V 5+ 
and Cr 6+ -based compounds would have even smaller 
band gaps than those previously synthesized, allowing for 
strong absorption in the solar spectrum. Indeed, there 
exist V 5+ and Cr 6+ non-perovskite oxides with optical 
gaps in the visible region, such as V2O5 (2.3 eV) |30j and 
SrCrC>4 (2.44 eV). [31| The discovery of more materials, 
especially in the highly tunable perovskite class, could 
be extremely valuable. It should be noted, however, that 
toxicity is a possible drawback of the widespread use of 
V 5+ and Cr 6+ -based compounds. |32| 

In what follows, we explore the band gaps of twelve hy- 
pothetical V 5+ and Cr 6+ double perovskites with rock- 
salt B-site ordering (Fig. GQi)— (Sr,Ba) 2 (Al,Sc,Ga)V0 6 
and (Sr,Ba)2(Mg,Ca,Zn)Cr06- From our calculations 
of phonon frequencies, elastic constants, and cation 
ordering energies summarized in the Appendix, these 
twelve hypothetical cubic double perovskites are likely 
metastablc at room temperature. Full analysis of the 
global stability of these structures would require double 
perovskite energies to be compared to those of numerous 
other phases (e.g., binary and ternary oxides). However, 
their metastability suggests that these compounds could 
in principle be synthesized, if not in bulk then via epi- 
taxial growth. 

In addition to being metastable, the predicted elec- 
tronic structure of these new V 5+ and Cr 6+ -based com- 
pounds (Fig. [TJd) confirms the intuition that smaller, 
more highly charged metal cations result in lower-energy 
conduction bands and correspondingly smaller gaps. 
However, the LDA gaps implied by Fig. [TJd are expected 
to be underestimated and therefore not quantitative, 



Corrections of DFT band gaps to experiment 

In the remainder of this article, we develop and apply 
an ad hoc correction to DFT and many-body perturba- 
tion theory gaps of d° A 2 B'B06 double perovskites, es- 
tablishing their optical gaps within ±0.20 eV as signifi- 
cantly lower than previous materials studied in this class. 
We fit the difference between computed and measured 
gaps of a set of d° perovskite-derived compounds mea- 
sured by Eng et al. [llj, using them as a test set for our 
ad hoc correction. The compounds measured by Eng et 
al. [lO ] include a wide range of elements, their electronic 
structures are qualitatively similar, and their optical gaps 
have been measured using a consistent procedure. Of the 
31 compounds measured, we compute 30 (shown in the 
top section of Table |TJ in LDA and G W @LDA, exclud- 
ing only NaNbC>3 because its large 40-atom unit cell [IH 
is currently a computational challenge for Go Wo- Be- 
cause of the additional computational expense associated 
with optimizing crystal structures with hybrid function- 
als, we compute only the twelve cubic compounds with 
HSE06. Structures are optimized in their experimentally 
determined room-temperature space groups with rock- 
salt B-site ordering. 

In Fig.m LDA (red), GoWo@LDA (blue), and HSE06 
(green) gaps are plotted versus the optical gaps ob- 
tained by Eng et al. [ll| Double perovskite compounds 
are grouped according to the more highly charged B- 
site cation, whose d states are primary contributors to 
the conduction band edge. Optical gaps are underesti- 
mated by LDA, and significantly overestimated by one- 
shot Go Wo OLD A and HSE06. As expected, compounds 
with smaller, more highly charged cations (e.g., Mo 6+ ) 
tend to have smaller gaps than those with larger, less 
highly charged cations (e.g., Ta 5+ ). An unsatisfactory 
feature of Fig.[2]is that for all B-site cations, deviations in 
computed gaps from experiment are large and vary from 
compound to compound. We find mean average errors of 
-0.74 ± 0.30, 0.92 ± 0.36, and 0.57 ± 0.25 eV for LDA, 
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FIG. 1: a) The double perovskite structure, with rock-salt B-site ordering, b) LDA densities of states (DOS) of SrTiO;?, 
Sr2ScVC>6, and Sr2CaCr06, in which the valence band maxima are defined as zero energy, c) LDA (light red), GoWoCSLDA 
(light blue), and HSE06 (light green) band gaps of twelve hypothetical V 5+ and Cr 6+ double perovskites with rock-salt B-site 
ordering. Estimated optical gaps (red, blue, green) based on the assumption that computational errors are equal to those of 
SrTi03. Estimated optical gaps (black) based on the modification of LDA gaps using linear corrections in and Ax (see 
text). Error bars (±0.20 eV) indicate uncertainty, as determined by the standard deviations of the fits for known compounds. 
The dashed line represents the optical gap of SrTiC>3. 



TABLE I: (Top) Perovskite- derived compounds with d° B-site metal cations considered in this work, with optical gaps taken 
from the experiments of Eng et al. [ll|] (Bottom) Newly-predicted V 5+ and Cr 6+ double perovskites. All estimated gaps are 
based on Eqn. [T] and have uncertainties of ±0.2 eV. 
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GqWo@LDA, and HSE06, respectively. Applying cor- 
rections to computed gaps of V 5+ and Cr 6+ compounds 
based on these mean average errors results in estimates 
consistent with our earlier estimates based on SrTiC>3. 
However, as this approach assumes that computed gap 
errors are randomly distributed, a clearer understanding 
of the systematic variation of the gap errors within the 
test set is needed. 

Past work on traditional semiconductors and insulators 
has shown that underestimates of experimental gaps by 
LDA tend to roughly scale as — , where eoo is the elec- 
tronic component of the static dielectric constant. [44[ 
That is, LDA tends to underestimate experimental gaps 
more severely as — increases. In contrast, for d° per- 
ovskites and double perovskites, Fig. [3^ shows that the 
LDA gap is closer to experiment as — increases. The 
fact that this trend in Fig. [3^ runs opposite to established 
intuition suggests that there is another important factor 
determining band gap errors in this class of compounds. 
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FIG. 2: LDA (red), G W @LDA (blue), and HSE06 (green) 
electronic gaps versus experimental optical gaps of com- 
pounds synthesized and measured in reference Com- 
pounds are grouped according to the more highly charged B- 
site cation, the primary component of the conduction band. 
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FIG. 3: a) Differences between optical and LDA gaps (Alda) versus ^- for 30 compounds synthesized and measured in 

reference [TJ]. b) Alda versus Ax (see text), c) Alda versus both — and Ax. Dashed lines and planes are least-squares fits 
to the data shown. Compounds are grouped according to the more highly charged B-site cation, the primary component of the 
conduction band. 



Band edge states of d° perovskite-derived compounds 
are known to vary in their degree of localization, particu- 
larly the conduction band minimum (CBM). While cubic 
ABO3 pcrovskites have CBM states consisting purely of 
B-site metal cation d character, double perovskitc CBM 
states also include varying amounts of oxygen 2p charac- 
ter, [ll| leading to a range of localization depending on 
fractional d and 2p content. In Figure|3jD, we demonstrate 
that these differences in band edge localization are par- 
tially responsible for the errors in computed band gaps. 
As a proxy for band edge localization, we compute the 
difference between the band edge Fock exchange energies 
(Ax) from the Kohn-Sham LDA wavefunctions associ- 
ated with the CBM and valence band maximum (VBM). 
For these compounds, because their VBM exchange en- 
ergies are all quite similar, a larger value of Ax indicates 
that the CBM is derived from more localized d states. 
Not surprisingly, the data show that a more localized 
CBM generally results in a more severely underestimated 
LDA gap. 

In light of Fig.Eb, the counterintuitive trend in Fig. [3^, 
can be rationalized. In this class of d° perovskite-derived 
compounds, those which have the smallest LDA gaps 
(and therefore the smallest values of ^-) are those with 
the most localized CBM states (and therefore the most 
severe underestimates of exchange). 

In Fig. [3};, we show the result of fitting the LDA gap 
errors in both j— and Ax simultaneously, obtaining the 
following expression to predict the optical gaps of d° 
perovskite-derived compounds: 



£ s (opt) = £ g (LDA)+0.148A x - 



8.62eV 



-+0.33eV±0.20eV. 

(1) 

Compared to our earlier corrections to computed gaps, 
this fit reduces the uncertainty of optical gap predictions 
significantly, to ±0.20 eV. Estimates based on Eqn. [T]of 
band gaps within the test set are shown in the top section 
of Table |U 



Predicted band gaps of new V 5+ and Cr 6+ -based 
compounds 

In Fig. nj. (black) and the bottom section of Table |H 
we apply these refined band gap estimates to the twelve 
hypothetical V 5+ and Cr 6+ double perovskites discussed 
earlier. Error bars indicate the uncertainty of the esti- 
mates (±0.20 eV), as determined by the standard devi- 
ations of the fits for known compounds. The results in 
Fig. [TJ; indicate that V 5+ and Cr 6+ double perovskites 
with rock-salt B-site ordering would likely have optical 
gaps spanning significant parts of the solar spectrum 
(1.7-2.4 eV for V 5+ and 1.1-1.9 eV for Cr 6+ ). While 
band gaps in this region do not guarantee strong solar 
absorption, d° perovskite-derived compounds are likely 
to absorb strongly because their valence and conduction 
bands consist primarily of oxygen 2p and metal d states, 
respectively. 

We note that our band gap estimates from this linear 
fitting procedure hold for known d° compounds outside 
the test set. [56j For cubic perovskite BaZrC>3, we predict 
a gap of 4.1 ± 0.2 eV, compared to the measured gap of 
4.0 eV. For ferroelectric perovskite LiNbOs, we pre- 
dict a gap of 4.2 ±0.2 eV, compared to the measured 3.78 
eV. [46( The difference between theory and experiment 
for LiNbOs is not surprising given that zero-temperature 
calculations tend to exaggerate the amplitude of ferro- 
electric distortions. [22j It is likely that our approach can 
be used to predict the optical gaps of not only d° per- 
ovskites, but of d° transition metal oxides in general. For 
rutile Ti0 2 , we predict a gap of 3.1 ± 0.2 eV, compared 
to the measured value of 3.05 eV. 
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The fact that GoWo@LDA overestimates the optical 
gaps of our test set of compounds (by an average of 
0.92 eV) may stem from several sources. First, GW 
neglects electron-hole interactions, which are expected 
to be non-negligible in these systems (on the order of 
0.1-0.2 eV in recent calculations of Ti02 [Hj]). Sec- 
ond, electron-phonon interactions, also neglected in this 
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work, can renormalize the optical gap. [19|, |48[ Third, 
it is possible that stronger correlations, inaccessible via 
the GW approximation, are quantitatively significant in 
these systems. 4|| 5^] Finally, a primary assumption of 
one-shot GoIUoCJLDA is that the LDA wavefunctions are 
a good approximation to the quasiparticle wavefunctions. 
Further studies, employing a self-consistent approach to 
GW, and/or more optimal starting points, may improve 
agreement with experiment. [HI We plan to explore these 
ideas in future work. 



CONCLUSIONS 

In summary, we have predicted the band gaps and en- 
ergetic stability of a class of as-yet unmade A2B'B06 
double perovskites containing B-site V 5+ or Cr 6+ . 
Through a combination of LDA, GoWo@LDA, and 
HSE06 calculations, and analysis of the experimentally 
determined optical gaps of related d° perovskitc-derived 
compounds, we have determined that a number of these 
hypothetical V 5+ and Cr 6+ compounds are likely to be 
metastablc at room temperature with respect to distor- 
tion and changes in B-site ordering, and to have optical 
gaps spanning significant portions of the solar spectrum. 
If compounds such as these could be synthesized, whether 
in bulk or cpitaxially, they would likely be of interest in 
light harvesting and energy conversion applications. 
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FIG. 4: a) The most unstable phonon frequencies at the zone 
center of a 2 x 2 x 2 perovskite unit cell for twelve hypothet- 
ical V 5+ and Cr 6+ double perovskites with rock-salt B-site 
ordering, with phonon modes in which B-site cations move 
( "translational" ) shown separately from others ( "rotational" ) . 
b) Computed phonon dispersion relations of Ba2ScVOs (left) 
and Ba2CaCrOe (right) with undistorted rock-salt B-site or- 
dering. Real frequencies indicate stability with respect to 
distortion. 
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APPENDIX: METASTABILITY OF PREDICTED 
COMPOUNDS 

In what follows, we comment on the stability of V 5+ 
and Cr 6+ double perovskites studied in this work. For the 
rock-salt B-site ordering of our twelve hypothetical V 5+ 
and Cr 6+ compounds, structures are relaxed in DFT- 
LDA within FrnZm space group symmetry. All elas- 
tic constants (i.e., matrix elements of the stiffness ten- 
sor) are found to be non-negative. The most unstable 
phonon frequencies at the zone center of a 2 x 2 x 2 



perovskite unit cell are shown in Fig. Hk, separated be- 
tween phonon modes that are associated with the move- 
ment of B-site atoms ( "translational" ) and those that are 
not ( "rotational" ) . Seven of the twelve compounds have 
only real phonon frequencies, indicating stability with 
respect to structural distortions. The five compounds 
with imaginary zone-center phonon frequencies exhibit 
soft octahedral rotation modes of a magnitude similar 
to SrTiOa (91i cm -1 in LDA), where these modes are 
suppressed well below room temperature. [5^ 



53j Two 



compounds (S^CaCrOg and S^ZnCrOg) have unstable 
distortions that break inversion symmetry; we note that 
ferroelectric oxides with small band gaps may be of po- 
tential interest for their ability to separate electrons and 
holes. [12|, |54| Fig. |4)d shows full phonon dispersion re- 
lations for Ba2ScV06 and Ba2CaCrO@, which indicate 
stability throughout the Brillouin zone. Due to their 
structural and chemical similarity, we expect qualita- 
tively similar behavior for the other hypothetical com- 
pounds. These results suggest that most of these hy- 
pothetical cubic double perovskites with rock-salt B-site 
ordering are likely metastable at room temperature. 
To assess whether rock-salt ordering is indeed the most 
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FIG. 5: a) The six distinct B-site orderings of A2B'BC>6 dou- 
ble perovskites in a 2 x 2 x 2 perovskite unit cell, b) Total ener- 
gies per formula unit, relative to the rock-salt B-site ordering 
of each compound and in the absence of symmetry-breaking 
distortions, c) The most unstable phonon frequencies at the 
zone center of a 2 x 2 x 2 perovskite unit cell for the distinct 
B-site orderings of Ba2ScV06 and Ba2CaCrOe, with phonon 
modes in which B-site cations move ( "translational" ) shown 
separately from others ("rotational"). 



stable arrangement of B-site atoms, the total energy 
associated with this ordering was compared to those 
of the five other distinct B-site orderings possible in a 
2x2x2 perovskite unit cell (Fig. [5^). In the absence 
of symmetry-breaking distortions, rock-salt ordering has 
the lowest total energy for all twelve combinations of el- 
ements (Fig. [S})). 

However, there is still the possibility of distortions 
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FIG. 6: The DFT-LDA total energy per formula unit of the 
ferroelectric (FE) columnar ordering of twelve hypothetical 
V 5+ and Cr 6+ double perovskites relative to rock-salt ordering 
(blue), and DFT-LDA gaps of both geometries (green, red). 



(translational or rotational) changing the relative ener- 
getics of these orderings. To explore this possibility, 
Ba 2 ScVC>6 and Ba 2 CaCrC>6 were subjected to more thor- 
ough stability analysis. Their most unstable phonon fre- 
quencies at the zone center of a 2 x 2 x 2 perovskite 
unit cell are shown in Fig. [5J; for all six B-site order- 
ings in the absence of distortion. In all cases, elastic 
constants are non-negative. Based on these results, we 
expect octahedral rotations to have little effect on the 
relative energetics, as they are equally present in all or- 
derings (and possibly suppressed at room temperature as 
in SrTiOs). Structures were then relaxed along all un- 
stable translational phonon modes along high-symmetry 
axes. For both Ba 2 ScVC>6 and Ba 2 CaCr06, the only 
structures whose energies compete favorably with those 
of rock-salt ordering are distorted variants of columnar 
ordering (Fig. [5^, top middle). In both cases, the struc- 
ture with the lowest total energy is the columnar struc- 
ture with ferroelectric polarization along the direction of 
the columns. Both compounds in this geometry have 
non-negative elastic constants and real phonon frequen- 
cies, meaning they are not expected to distort further. 
In Fig. [6j the rock-salt and polarized columnar B-site 
orderings of all twelve hypothetical compounds are com- 
pared in both total energy (blue) and DFT-LDA band 
gap (green, red). For three of the twelve compounds 
(Sr 2 AlV0 6 , Ba 2 AlV0 6 , and Sr 2 MgCr0 6 ), rock-salt B- 
site ordering is the most stable structure. For the re- 
maining nine, a few points arc worth noting. First, it 
is possible that internal distortions would be sufficiently 
suppressed at room temperature as to favor rock-salt or- 
dering. Second, the computed DFT-LDA gaps of the two 
competing structures are quite similar in most cases, sug- 
gesting that even if present, this change in B-site ordering 
and distortion would do little to change measured band 
gaps. Third, even if it is not the lowest-energy structure, 
rock-salt B-site ordering could in principle be achieved 
via epitaxial growth in the [111] direction. 
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